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ABSTRACT 

We study the specific star formation rate (SSFR) and its evolution at z >, 4, in mod- 
els of galaxy formation, where the star formation is driven by cold accretion flows. 
We show that constant star formation and feedback efficiencies cannot reproduce the 
observed trend of SSFR with stellar mass and its observed lack of evolution at 2; > 4. 
Model galaxies with log(M*) ^ 9.5 Mq show systematically lower specific star forma- 
tion rates by orders of magnitudes, while massive galaxies with M* ^ 5 x 10^*^ M© 
have up to a n ord er of magnitude larger SSFRs, compared to recent observations by 
IStark et"al] (|2009t ). To recover these observations we apply an empirical star forma- 
tion efficiency in galaxies that scales with the host halo velocity dispersion as ex 1/a^ 
during galaxy mergers. We find that this modification needs to be of stochastic na- 
ture to reproduce the observations, i.e. only applied during mergers and not during 
accretion driven star formation phases. Our choice of star formation efficiency during 
mergers allows us to capture both, the boost in star formation at low masses and the 
quenching at high masses, and at the same time produce a constant SSFR-stellar mass 
relation at z ^ 4 under the assumption that most of the observed galaxies are in a 
merger triggered star formation phase. Our results suggest that observed high-z low 
mass galaxies with high SSFRs are likely to be frequently interacting systems, which 
experienced bursts in their star formation rate and efficiency (mode 1), in contrast 
to low redshift z ^ 3 galaxies which are cold accretion-regulated star forming systems 
with lower star formation efficiencies (mode 2). 

Key words: galaxies: formation - galaxies: high-redshift - galaxies: interactions - 
galaxies: starburst 



1 INTRODUCTION 

The cosmic star formation rate in the universe is shaped 
by the physical conditions that govern the growth of galax- 
ies. An initial, probably short-lived, t ~ 10^ yr, episode of 
population Ill-d ominated star formation at z ~ 15 (e.g. 
iMaio et al.l [2010|) is followed by a steady increase in cos- 
mic population II/I star formation activity that re aches 
a peak around z ~ 2 (| Hopkins fc BeacomI 120061 ) ■ co- 
inciding with a peak in the merger activity of galax- 
ies (e.g. iKhochfar fc BurkertI I2OOII : IConselice et"ai] I2OO8I : 
iKhochfar fc Silkll2009ah . The steady decline at lower red- 
shifts is quite strong, however, allowing for ~ 50% of the 
presen t-day st ellar ma ss to have be en formed between z = 
and 1 jPerez-Gonzalez et al.ll2008l '). 

Within the hierarchical ACDM framework, the initial 
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increase in the cosmic star formation is closely linked to 
the growth of the hosting dark matter haloes. The short 
cooling time-scale of halo gas favours efficient fuelling of star 
formation in galaxies, in particular at high redshi fts and 



in dark matter haloes with Myjr ^ 10 Mq (e.g. iBinnev 



1977]: [Sm3 I1977I : iRees fc Ostrlkeil Il977l : iBirnboim fc Dekel 
2003'). More massive haloes at z > 2 however, reveal efficient 
feeding of their central galaxies as well. Cold streams of gas 
that enter the virial radius along cosmic filaments penetrate 
through the diffuse h ot halo gas and reach the potentia l 
minimum of the h alo JKeres et al.ll2005l : [Ocvirk et al.ll2008l : 
iDekel et al.ll2009l ). Massive galaxies at z ~ 2 show gas-to- 
stellar mass ratios of order unity suggesting that cold flows 
could in deed be a very effic i ent means of provid ing gas to 
galaxies (|Daddi et al.ll2010al : iTacconi et akfeoiOl ). 

At z < 2, cold streams in massive haloes cease to 
exist and feedback effects co ntrol the cosmic star forma- 
tion (e.g. ISchave et al.l I2OI0I . see however Bouche et al. 
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2009). Different feedback sources are also held responsi- 
ble for shaping the mass function (e. g. Dekel fc Silk 198(: 



Benson ct al. 2003:'Cr oton et al.ll200alKhochfar et al.ll2007 : 
Oppenheimer ct al. 2010|') or the colou r s of galaxies (e.g. 
Kang et al.1 l2005l: ICattaneo et al.1 l2006l : iDekel fc BirnboimI 



20061 : iKhochfar fc Ostrikeijl2008l ) 



The importance of cold accretion for the star forma- 
tion rate of galaxies has recently been the subject of vari- 
ous semi-analytic studies, which show that the evolution of 
the specific star for mation rate (SSFR) of ga,laxies can be 
reproduced at 2 < 4 jKhochfar fc Silkll20q9bi : iBouche et al.1 
I2OO9I : button et alll201(]l : iKang et al.ll2010l 'l. However, these 
model s fail to reproduce the observed SSFR of galaxies at 
z > A l|Bouche et al.ll2009r). where in partic u lar the SSFR is 
appro ximately constant (jStark et al.l l2009l : iGonzalez et al.l 
l2010f ). One option is that a second mode of star forma- 
tion enters at high redshift, perhaps driven by mergers 
l|Genzel et al.ll201Cll : iDaddi et aklbOlObl ). in contrast to the 
canonical cold gaseous disk instability model that elegantly 
reproduces the Schmidt-Kennicutt star formation relation 
at z < 2. 

In this Letter, we address the question whether it is pos- 
sible to reconcile observations of specific star formation rates 
at 2 > 4 with simple models of accretion-driven star forma- 
tion at high-z. We focus on the question how star formation 
needs to be modulated to accomplish this task and derive 
empirical scaling relations for the star formation efficiency 
in galaxies at high z. 



2 MODEL 

The basic modelling approach we follow is similar to 
iNeisteui fc WeinmannI ||2010t ) and iBouche et al.1 l|2009l ). In 
this approach the evolution of a galaxy within a dark mat- 
ter halo is followed by calculating the accretion rate of cold 
gas, the associated star formation rate and the effect from 
supernovae feedback for a grid of dark matter haloes with 
varying masses starting at a fiducial redshift. We approx- 
imate the accretion rate of cold gas onto central galax- 
ies by Mace = fbMoM, with fb — 0.165 as the univer- 
sal baryon fractio n and Mdm the dark matter halo mass 
IjDekel et al.ll2009l ). To simplify matters even more, we ap- 
proximate the dark matter growth rate by an analytic fit 
Mdm (z) « 6.6/-^MiV^(l + 2)^'^^ Mq yr-i (|Neistein et al.1 
120061 ). Here M12 is the dark matter halo mass in units of 
lO^'^ M0. Once the gas accretion rate as a function of time is 
known, one can calculate the gas mass in the central galaxy 
of a given dark matter halo. We assume that the gas set- 
tles onto a rotationally-suppo rted disc hosted by an isother- 
mal halo using the model of IMo et al.l l|l99a) . and convert 
the gas into stars using a global Schmidt-Kennicutt law 
M* = OitMgas/tdyn (jKennicutt 1998). The fiduci al value for 
the s tar formation efficiency is set to a, — 0.02 (|Kennicuttl 
Il998l'l and the dynamica l time of the disc, tdy„ = O.lRyir/Vc 
( Kauffmann et al.lll999l l. Here Rvir is the halo virial radius 
and Vc its circular velocity. Besides models with constant 
a, we also investigate models with more general efficiencies 
that depend on the halo velocity dispersion. The feedback 
effects from supernovae onto the cold gas rese r voir a re in- 
cluded using the approach in iKauffmann et aU l| 19991 ') with 
MsN = eEsNtniMFMt/Vc- The energy per supernova is 



EsN = 10^^ ergs, miMF = 1/150 Mq^ is the number of 
type-II supernovae per solar mass stars formed assuming 
a Chabrier IMF, M* the star formation rate and e is the 
feedback efficiency for which we assume a fiducial value of 
e = 0.2. The gas reheated by supernovae feedback con- 
tributes to the hot halo and we assume that it does not 
cool down by 2 = 4. Numerical simulations suggest that 
most gas i n galaxies indeed never reheats before joining it 
(JKeres et al. 2009 ). In the following we show our results for 
solving the differential equations for the evolution of the stel- 
lar mass M, , the cold gas mass Mgas = Mace — M^ — Msn, 
and star formation rate M* , by explicitly following the dark 
matter growth Mdm as a function of time using the fitting 
formulae for the average growth rate given above. We start 
at 2 = 20 with a grid of dark matter haloes ranging from 
10^ — 10^" Mq and set the initial conditions of the galaxies 
in the haloes to M, = Mgas = 0. We show in the follow- 
ing the properties of individual central galaxies that live in 
dark matter haloes with such average growth histories. Thus 
the galaxy properties refer to galaxies with average growth 
histories. In the following we assume a standard cosmology 
with {Q.^,Q.K,h,Gi.) = (0.3,0.7,0.7,0.8). 



3 RESULTS 

In Fig. [T] we show the predicted SSFR of galaxies for the 
fid ucial model and co mpare them to the observational data 
bv lStark et"ai] l|2009D . 

Our resul ts show the same tre nd as the non-accretion 
floor model of IBouche et al.l ((20091). Applying an accretion 
floor of Mdm = 10^^ M©, below which no cold accretion 
occurs, gives good agreement with observed data at 2 < 3, 
but suppresses the formation of enough low mass galaxies 
at 2 > 5. For this reason we do not apply an accretion floor, 
and allow gas accretion in halos of all mass. The trend in 
SSFR is in stark contrast to the observed data. Most strik- 
ingly, low mass galaxies of M» ;$ 5 x 10® M© show orders 
of magnitudes too low star formation rates, while massive 
galaxies with M» ^ 5 x 10^*^ M© show an order of magni- 
tude too high star formation rates. Independently of these 
failures, the time-scales of the inverse SSFRs at 2 = 4 are 
smaller than the age of the universe, indicating efficient early 
formation. Fig. [T] also highlights that the model SSFR-Af« 
relation is significantly evolving from 2 = 4 to 2 = 6, in- 
dicated by open and filled circles, respectively. In contrast 
the observations shown in the same figure indicate almost 
no evolution. 

Before we continue comparing model results to obser- 
vations it is worth pointing out potential observational bi- 
ases. The main source of uncertainty are the observed star 
formation rates. In Fig. [T] w e show SSFRs th at are not 
corrected for dust extinction IjStark et al.ll2009l ). Including 
dust extinction results in higher star formation rates, shift- 
ing the SSFR-M* relation to larger values in SSFR, thus 
worsening the discrepancy between models and observa- 
tions in particular at the low mass end, if these galaxies 
were to host significant amounts of dust. Another important 
source of unc ertainty is the com pleteness of the observed 
samples. The IStark et al.l l|2009l ) sample is co mplete above 
log(gF fl) = (0.25,0.41,0.54) for 2 = (4. 5. 6). [Stark et all 
(J2009) argue that at low stellar masses their sample is not 
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Figure 1. The specific star formation rate as a function of stel- 
lar mass. Each individual line shows the growth history of one 
central galaxy within a dark matter halo under the assumptions 
of the fiducial model. The initial halo masses at z = 20 span 
a range from 10^ — 10^*^ Mq. The solid black circles and open 
circles show the model SSFR at z = 4 and z = 6, respectively. 
The red triangles and blue squares show the observed data from 
IStark et al.l l|2009l ) at these redshifts, uncorrected for dust extinc- 
tion and assuming an exponential decaying star formation rate 
with T = 100 Myr. Typical error s are 0.3 dex. Light b l ue tr i- 
angles show the z = 7 data from ISchaerer fc de Barro j 1 120101) . 
Dashed lines show lines of constant star formation rate for 500, 
100, 10, 3.5 and 1.8 M p/yr. Latter two co rrespond to the reported 
completen ess limit in [ S tark et al.| 1120 08) at z = 6 and 4, respec- 
tively. The ISchaerer fc de Barroj 1(20101) sample is a combination 
of three different samples and we refer the reader to their paper 
on a discussion on the observational limits and incompleteness of 
their sample selection. The solid horizontal line marks the inverse 
age of the universe at z = 4. 



complete, and misses potentially low star forming galaxies. 
If indeed a large population of galaxies with star forma- 
tion rates below these limits exists, the average derived SS- 
FRs are over-estimates. However, independent of such cor- 
rections, the upper envelope of the SSFR-M, can be trusted, 
and poses particularly at low masses a challenge to models. 
We tested the behaviour of the SSFR for various star 
formation and feedback efficiencies. In general, increasing 
a, leads to an increase in the stellar mass of gala:xies within 
halos, and a lower SSFR. The latter is directly related to the 
scaling of the accretion rate with halo mass Mace oc M}^]^j 
and the lower reservoir of remaining gas in galaxies. The 
same trend is found for reducing the feedback efficiency e. 
SSFRs drop as galaxies form more easily and now live in less 
massive halos. Our results suggest that a constant boost in 
star formation efficiency is not going to improve the low 
SSFR for galaxies with M , < 10^ M© . What is needed is a 
boost in a* that is not constant, but of a stochastic nature. 
Within the CDM framework, galaxy mergers are a natural 
candidate to provide such a behaviour. To recover the strong 



Figure 2. Same as Fig. [T]for a model with a, (a) = olq (o'o/cr)'' 
during mergers, with ctq = 85 km/s, 7 = 3 and ao = 0.02. During 
non-merger episodes o = oq- We here assume that all galaxies 
undergo mergers at z; = 4, 5, 6, 7, 8, 9, 10. In reality galaxies will 
experience mergers at various different redshifts which will result 
in a scatter in the SSFRs. Thus the relation we are showing here 
represents the maximum SSFR-M* relation for individual galax- 
ies. Galaxies not going through mergers will be lying closer to the 
evolutionary tracks shown in Fig. [T] The z = 6 model-relation 
looks virtually identical to the one at z = A, which is why we 
omit showing it in the figure. 



observed mass-dependence of the SSFR, we propose a rela- 
tion in which the star formation efficiency during mergers 
is a power-law of the halo velocity dispersion. We adopt the 
following relation during mergers: 



Q. = ao 



m^ 



(1) 



with ao = 0.02 as the fiducial value observed locally and ao 
a characteristic velocity dispersion. The values of 7 and ctq 
are treated as free parameters to obtain the observed scal- 
ing in the relation between the SSFR and the stellar mass 
at_z> 4. Our best-fit values matching the z = 4 lStark et al.l 
(J2009|) data are ctq ~ 85 km/s and 7 = 3. A si milar scaling 
has been recently suggested in lSilk fc NormanI (1200.9) based 
on a porosity model for the inter-stellar medium. Within 
their model black hole feedback triggered by mergers plays 
a prominent role. The authors argue that black holes in 
low mass galaxies are not massive enough to shut down 
star formation during their init ial growth phase, but in- 
stead can possible help boo st it (|Antonuccio-Delogu fc SilH 
I2OO8I : ISilk fc NormanI [2009|) , whereas massive black holes 
are able to prohibit sta.r form ation via AGN-feedback (e.g. 
ISilk fc Rees|]l998l : lKinell2005l ) . 

We implement the stochastic nature of mergers in a very 
simplified manner to highlight their effects on the SSFR. At 
redshifts z = 10, 9, 8, 7, 6, 5 and 4, we assume that every 
galaxy has an instantaneous merger. We assume that the 
burst is lasting 50 Myr with an efficiency given by Eq. [l] 
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We furthermore neglect the stellar mass contribution of the 
merging partner to the remnant. Our model is clearly a sim- 
plification, and we tested the influence of our assumptions 
on the results by letting mergers occur randomly during 
the evolution of galaxies, and changing the duration of the 
merger episode as well as adding stellar mass from a merger 
partner. It turns out that the exact frequency of mergers is 
not important, as long as there is enough time for gas accre- 
tion to happen after the merger. Given the high accretion 
rates at z > 4 merger-triggered star formation can be fre- 
quent in our model, as long as the merger duration is not too 
long (< 100 Myr). Adding mass from merger partners only 
shifts the relation slightly, which can be compensated for by 
adjusting ctq- Outside merger episodes, we apply the fiducial 
star formation efficiency, i.e. «» = ao and follow the evolu- 
tion of our model galaxies. The evolution of the SSFR is pre- 
sented in Fig. [2] We find good agreemen t between the mod - 
elled SSFR- A'/* relation and the data of lStark et"ai1 l|2009l ). 
In addition we find that the our SSFR-M* relation does not 
change much with redshift, if we assume interacting galax- 
ies are dominating the sample. As we pointed out above, 
the observed sample might be missing a population of low 
star forming galaxies. Such galaxies would be non-merging 
galaxies within our model, which show low SSFRs closer to 
the fiducial model shown in Fig. [l] To highlight the non- 
evolution of the SSFR we compare in Fig. [3] our model pre- 
dictions for individual galaxy masses with the observations 
in the literature. We fi nd good agreement w ith the data of 
IStark et all (|2009l ') and lGonzalez et all (|2010l ') if we consider 
model galaxies between 10^ and 5 x 10® Mp). However, the 
data of 'Schaerer fc de Barroj (|2010l ) and lYabe et"ai] (|2009l ') 
show larger values in SSFR. Most of the discrepancy be- 
tween these and the former data sets is due to different star 
formation rate estimates that the authors assume, reflect- 
ing e.g. their different choice of dust extinction corrections 
and star formation histories. Interestingly, the latter sur- 
veys also predict no or only a week evolution from z = 4 
to z = 7. Within our model we can achieve higher SSFRs 
to be in agreement with this data by adjusting our free pa- 
rameter ap. Our r esults suggest, that if indeed the sample of 
I St ark et al.l (|2009|) is complete in terms of SSFR sensitivity, 
most if not all low mass high star forming galaxies have to 
be merger induced. 



4 DISCUSSION AND CONCLUSION 

In this Letter, we present a study of the SSFR in ACDM 
models in which galactic star formation is driven by cold gas 
accretion onto a central galaxy, that parallels the growth of 
the hosting dark matter halo. Such models have been shown 
to reproduce the observed star formation rates of galaxies 
successfully at z ;S 3 (e.g. Bouchc ct al. 2009), but fail at 
high redshifts (see however iLacev et al.l [2010|) . The rise in 
SSFR at 2 ^ 2 is due to the increased gas fraction. However, 
a new ingredient is required at higher z. The two observa- 
tional constraints of a high SSFR for low mass galaxies and 
non-evolution of the SSFR-stellar mass relation at z > 4 
cannot be met by a simple constant increase in the star for- 
mation efficiency a* at z > 4. However, an occasional boost 
in a, triggered by mergers in between periods of low star 
formation efficiency can provide high SSFR. Galaxy mergers 
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Figure 3. The evolution of the SSFR for the merger trig- 
gered boosts in at . For details on the low redshift data see 
iBouche et"all ||2009|) At z < 7, typical errors in SSFRs are ~ 0.3 
dex l| Gonzalez et al.l |2010|) , indicated by th e error-bar on the 
top left. We show the average data for the IStark et al.l 1120091 ) 
sample limiting ourselves to galaxies with M» ^ 10® Mq and 
star formation rates above the cited completeness limit. The 
iGonzalez et al.l |2010) data is the average for their sample, which 
has a medium m ass of A/« ~5 X 10® Mq. The lower and up- 
per ISchaerer &: de Barrog l|2010l ) points refer to ga laxies with 
log(M «) ~ 9.5 and log(M«) ~ 8, respectively. The lYabe et al.l 
||2009|) data includes galaxies with masses from 10* to 10^^ Ts/Iq 
with a median value of 4.1 X 10® Mq. We show model results for 
SSFRs averaged over mass bins as lines. We here assume that the 
haloes within a mass bin contribute equally. 



are known to trigger starbursts even in minor m erger inter- 
actions {M1/M2 < 10, Ml ^ Ah) (e.g. IPeiranTet al. 2010 ). 

Observationally, merger-induced star formation at low 
redshifts 2: < 2 i s not the main driver of star formation 
(|jogee et al.ll2009l ). This is consistent with the growth of 
dark matter haloes in simulations, which takes place through 
signi ficant smooth acc retion events and not major mergers 
(e.g. iGenel et al.ll2010l ). Our results indicate that the situ- 
ation at 2 > 4 is rather different. Merging time-scale esti- 
mates using dynamical friction predict to first order a depen- 
dence on the mass ratios of the mergin g partners only, once 
expressed in units of t he Hubble time |Springel et al.l I2OO1I : 
iBovlan-Kolchin et alj I2OO8I : I Jiang et al.l I2OO8I '). Aii equal 
mass merger thus will take place approximately a factor of a 
few longer at 2: = 2 than at 2 = 6. This behaviour supports 
frequent interactions at 2 > 4 that can trigger star forma- 
tion while at lower redshifts star formation takes mainly 
place during quiescent periods in between mergers. 

While mergers certainly can boost star formation, 
AGN-regulated starbursts in mergers could provide a nat- 
ural explanation for the non-evolution of the SSFR-stellar 
mass relation at 2 > 4. If the effects of AGN-feedback are r e- 
lated to the dark matter potential (jBooth fc Schavell2010f ) a 
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scaling of the star formation efficiency with the halo velocity 
dispersion can effectively produce a mildly or non-evolving 
relation. We show that this is indeed the case for the empir- 
ical scaling that we adopt in this paper with a* oc a~^ dur- 
ing mergers. A more physical model might appeal to AGN 
growth and feedback which are much less efficient in low 
mass galaxies, thereby inhibiting quenching, while in mas- 
sive galaxies, the central black holes are already fully grown 
and are responsible for quenching the SFR. 

Enhanced merging at high redshift provides the SSFR 
boost that also lends itself to generating enhanced a/Fe in 
the resulting massive galaxies. We predict that there should 
be rare low mass a/Fe "refugees", perhaps companions of 
massive ETGs, that have avoided the final merging fate and 
gas blow-out, but nonetheless carry chemical traces of their 
enhanced SSFR history. 

The role of AGN remains to be elucidated. Quenching 
of star formation is commonly attributed to AGN. This may 
be the case for massive galaxies at high redshift. However 
AGN may also play a role in boosting star formation in the 
low mass systems where the SSFR is enhanced. A future test 
of the role of AGN will be to examine the residuals of a sam- 
ple with measured AGN accretion rates and star formation 
rates in order to see whether for example the Eddington ra- 
tio correlates (boosting) or anti-correlates (quenching) with 
SSFR. 

Our results suggest that the majority of galaxies with 
M* ^ 10^ M0 at 2 > 4 are interacting systems. This pre- 
diction can be tested with future upcoming observational 
missions, and should provide a strong test on the early 
build-up of galaxies. In a follow-up paper, we will inves- 
tigate the individual contribution from merger triggered 
star-bursts and AGN with respect to accretion-driven star 
formation in the context of a full semi-analytic model. 

The authors would like to thank Dan Stark for kindly 
providing the observational data in electronic form and use- 
ful comments on the draft. SK acknowledges support from 
the the Royal Society Joint Projects Grant JP0869822. 
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